"Doubly- magic" conditions in "magic" -wavelength trapping of ultracold alkalis 
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In experiments with trapped atoms, atomic energy levels are shifted by the trapping optical 
and magnetic fields. Regardless of this strong perturbation, precision spectroscopy may be still 
carried out using specially crafted, "magic" trapping fields. Finding these conditions for particularly 
valuable microwave clock transitions in alkalis has so far remained an open challenge. Here I 
demonstrate that the microwave clock transitions for alkalis may be indeed made impervious to 
both trapping laser intensity and fluctuations of magnetic fields. I consider driving multiphoton 
transitions between the clock levels and show that these "doubly-magic" conditions are realized 
at special values of trapping laser wavelengths and fixed values of relatively weak magnetic fields. 
This finding has implications for precision measurements and quantum information processing with 
qubits stored in hyperfine manifolds. 

PACS numbers: 37.10.Jk, 06.30.Ft 
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Compared to spectroscopic beam experiments, trap- 
ping cold atoms and molecules removes Doppler shifts 
and increases interrogation time thereby dramatically en- 
hancing spectral resolution. This improvement comes at 
a price: trapping optical fields strongly perturb atomic 
energy levels - transition frequencies are shifted away 
from their unperturbed values. In addition, the underly- 
ing Stark shift is proportional to the local intensity of the 
trapping lasers; the shift is non-uniform across the atomic 
ensemble and it is also sensitive to laser intensity fluctu- 
ations. So trapping seems to be both advantageous and 
detrimental for precision measurements. This dilemma 
is elegantly solved using so-called "magic" traps [lj]. At 
the "magic" trapping conditions two levels of interest are 
shifted by exactly same amount by the trapping fields; 
therefore the differential effect of trapping fields simply 
vanishes for that transition. 

The idea of such "magic" trapping has been crucial for 
establishing a new class of atomic clocks [2|, the optical 
lattice clocks. Here atoms are trapped in optical lattices 
formed by counter-propagating laser beams; the lasers 
operate at the magic wavelength. In these clocks one 
employs optical transitions in divalent atoms, such as Sr 
and Yb. Finding similar magic conditions for ubiquitous 
alkali-metal atoms employed in a majority of cold-atom 
experiments remains an open challenge. Especially valu- 
able are the microwave transitions in the ground-state hy- 
perfine manifold. Identifying magic conditions here, for 
example, would enable developing microMagic clocks |3| : 
lattice clocks operating in the microwave region of the 
spectrum. In addition, the hyperfine manifolds are used 
to store quantum information in a large fraction of quan- 
tum computing proposals with ultracold alkalis. Finding 
magic conditions would enable a decoherence-free trap- 
ping for this important realizations of qubits. 

The clock transitions in divalent atoms are between 
non-magnetic states; this removes sensitivity to magnetic 
fields. For alkalis, however, an additional piece of the 



puzzle is that the clock/qubit states are sensitive to both 
optical and magnetic fields. One needs to eliminate the 
sensitivity of transition frequency v to both perturba- 
tions simultaneously. This problem is solved here. We 
will require that the clock transition is insensitive to 
both Stark- and Zeeman-induced perturbations (we will 
use the "doubly-magic" qualifier for such trapping con- 
ditions). 
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FIG. 1: (Color online) Left panel: Zeeman effect for the 
hyperfine manifold in the ground state of / = 3/2 iso- 
topes of alkalis. Two clock levels \F' = 2,M' F = +1) and 
\F — 1,Mf = —1) are shown in black. Right panel illus- 
trates geometry of laser-atom interaction: degree of circular 
polarization, angle 0, and laser wavelength may be varied. 

First steps in identifying magic conditions for hyper- 
fine transitions in alkalis have been made in Refs. [4j-l6|. 
That work focused on eliminating sensitivity to laser in- 
tensity It, by tuning the laser frequency to its magic value 
oj m : SvUuyj) = regardless of the value of II ■ The pro- 
posals 0, Q have neglected the effect of magnetic fields 
and focused on B-field sensitive Mp ^ states. So while 
the trapping would be Stark-insensitive, the states would 
decohere due to coupling to stray B-fields. Apartial solu- 
tion to this problem was discussed in Refs. [JQ: moving 
to the Mf = states eliminates sensitivity to Zeeman 
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shifts to the leading order. Yet one needs to apply a 
bias magnetic field of a specific value making the condi- 
tions "magic" for a given trapping laser wavelength. As a 
result, the transitions remain Zeeman-sensitive through 
the second-order effects; numerical estimates show that, 
unfortunately, the residual B-ficld sensitivity would pre- 
clude designing a competitive clock. 

Clearly, the sensitivity to B-fields has to be addressed. 
We require that at the magic B-field dv/dB{B m ) = 0. 
Such conditions occur, for example, for a two-photon 
\F' = 2,M' F = +1) \F = 1,M F = -1) transition in 
87 Rb at the field of about 3 Gauss. The relevant Breit- 
Rabi diagram is shown in Fig. [TJ The two clock levels 
are highlighted: the existence of the magic B-field may 
be inferred visually. Experimentally, such magic condi- 
tions have been proven instrumental for performing col- 
lisional studies in a magnetically-trapped Bose-Einstein 
condensates of 87 Rb @ and studying decoherence of cold 
cloud of Rb near microchip 10]. Similar ideas are es- 
sential to realizing atomic clocks on a microchip ll|, [12 1 . 
My present work adds the optical fields to the mix: it 
turns out that the unified Stark-Zeeman description is 
non-trivial due to an interference of the two effects. 

Formalism — We are interested in the clock tran- 
sition between two hyperfine states \F' = I + 1/2, M' F ) 
and \F = I — 1/2, Mp) attached to the ground electronic 
nSi/2 state of an alkali-metal atom (/ is the nuclear spin). 
Here and below we denote the upper clock state as |F') 
and the lower state as \F) (see Fig. [TJ). 

We focus on the M' F = —Mp transitions. For these 
transitions, the electronic ^-factors of the two states are 
the same (see Fig. [TJ). Then the bulk of the Zeeman 
shift of the transition frequency goes away and the linear 
Zeeman effect is determined only by the nuclear (/-factor 
gi = 1/1 /i ri uc/ Hn, where fj,N is the nuclear magneton. 
The residual linear shift is compensated by the second- 
order Zeeman correction, quadratic in the B-field. This 
leads to the magic value of the B-ficld 
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where \i e is the operator of magnetic moment of electron. 
The second-order estimate, Eq.([TJ, is a good approxi- 
mation: for 87 Rb it gives 3.25 G, while the "all-order" 
Breit-Rabi analysis yields B m = 3.228917(3) G, Ref. 0]. 
Values of B m for Rb and Cs isotopes are tabulated in 
Table |U Notice that the fields for all tabulated transi- 
tions are relatively weak and can be well stabilized using 
existing technologies 

Fixing magnetic field at its magic value accomplishes 
the Zeeman-insensitivity of the clock transitions. Now 
we would like to additionally remove the Stark sensitiv- 
ity to intensity of trapping laser fields. We consider the 
following setup shown in Fig. [TJ An atom is illuminated 
by a laser light, with a certain degree of circular polar- 
ization A. At the same time, a bias magnetic field is 

TABLE I: Values of "magic" B-fields and "magic" wave- 
lengths. 
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s, Rb, / = 3/2, 


fdock = 6.83 GHz 


|2,1>- 


>|1, 


-1) 3.25 


806 nm° 






85 Rb, / = 5/2, 


^lock = 3.04 GHz 


[3, 1> - 


> % 


-1) 0.359 




|3,2> - 


* I 2 , 


-2) 1.15 


479-658; 797-878 






133 Cs, / = 7/2, 


f dock = 9.19 GHz 


|4,1>- 


> |3, 


-1) 1.41 




[4, 2) - 


♦ |3, 


-2) 3.51 


906-1067; 560-677 


|4,3>- 


* |3, 


-3) 9.04 


898-1591;863-880; 512- 



"nearly doubly-magic 



applied at an angle to the direction of laser propaga- 
tion. The B-field is fixed at its magic value. This is a 
basic building block for optical trapping. For example, 
an optical lattice (standing wave) may be formed by two 
counter-propagating lasers of the same wavelength. 

In a laser field, both clock levels are shifted due to 
the dynamic (i.e., laser- frequency-dependent) Stark ef- 
fect [13( • I derived the followin g fo rmula for the differen- 
tial shift of the clock frequency [18j, 
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Here El is the amplitude of laser E- field, II oc E\ . The 
quantities (3 s , f3 a , and a a are the scalar and vector (axial) 
polarizabilities (see below) . An important fact is that all 
these polarizabilities depend on the laser frequency, uj. 
By tuning the laser frequency we require that the com- 
bination in curly bracket becomes zero. At that "magic" 
point, the differential clock shift vanishes independently 



of the laser intensity: Sv c \ oc ^(uj m ) = 0. 

What is the difference between the polarizabilities a 
and j3? We are considering the Stark shift of hyperfine 
levels attached to the same electronic state. To the lead- 
ing order, the shift is determined by the properties of the 
underlying electronic state (polarizability a). However, 
because the electronic state for both hyperfine levels is 
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the same, the levels are shifted at the same rate[19| and 
we need to distinguish between the two hyperfine levels. 
An apparent difference between the two clock levels is 
caused by the hyperfine interaction (HFI), and the rigor- 
ous analysis involves so-called HFI-mediated polarizabil- 
ities, /?. Lengthy third-order (two dipole couplings to the 
laser field and one HFI) expressions for these polarizabil- 
ities may be found in Ref. [4|. 

Continuing with our discussion of the clock shift ©, 
I would like to stress an unconventional origin of the 
last contribution. Its full form for an arbitrary B-field 
is A cos 9 Mf'91 1 ^ a nS 1/2 x {B/B m ). The term arises 
due to an interference between Stark and Zeeman in- 
teractions. Qualitatively, the vector contribution to the 
Stark shift has the very same rotational properties as the 
Zeeman coupling (both are vector operators) . These op- 
erators, in particular, couple the two hyperfine manifolds. 
Consider shift of the \F', Mp) level. The Zeeman oper- 
ator couples it to the \F, Mf) intermediate state, and 
then the vector Stark shift operator brings it back to 
the \F',Mf) level thereby resulting in the energy shift. 
This cross-term is of the same order of magnitude as the 
other two terms in Eq. ([2]) and has to be included in the 
consideration. 

Now we can find magic wavelengths by numerically 
evaluating atomic polarizabilities entering Eq. ([2]). To 
this end I used a blend of relativistic many-body tech- 



niques of atomic structure, as described in [14]. To 
improve upon the accuracy, high-precision experimental 
data were used where available. To ensure the quality of 
the calculations, a comparison with the experimental lit- 
erature data on static Stark shifts of the clock transitions 
was made. It is one of the problems where a consistent 
treatment is important and less sophisticated estimates 
may fail even qualitatively, see [|| for a discussion. 

Results - I will present the results of the calcula- 
tions in the following form. Since the magic condition 
corresponds to Sv c \ oc ^(uj m ) = 0, we may recast Eq. ([2]) 
into 
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The r.h.s. of this equation depends on the laser frequency, 
while the l.h.s. does not. Moreover, |Acos(9| < 1, there- 
fore the magic conditions would exist only if for a given 
u) the r.h.s. is within the range — \M' F \ and \M' F \. 

Doubly-magic trapping of 133 Cs is analyzed in Fig. [5] 
This atom is metrologically important. For the past four 
decades the SI unit of time, the second, has been defined 
as a duration of a certain number of periods of radiation 
corresponding to the transition between the two hyper- 
fine levels of the ground state of the 133 Cs atom. Cs 
clocks serve as primary frequency standards worldwide 
and there is a substantial investment in the infrastructure 
supporting these clocks. The most accurate Cs clocks are 
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FIG. 2: (Color online) Magic conditions for 133 Cs. A de- 
pendence of the product M F i Acosd on trapping laser fre- 
quency (in atomic units) is plotted. The shaded regions are 
bound by — \M F >\ and +jM F /| lines. Magic trapping for a 
\F' = 4, M' F ) -¥ \F = 3, -M' F ) clock transition is only pos- 
sible when the computed curve lies inside the corresponding 
shaded region. 



fridge-sized fountain clocks (see, e.g., Ref. [l5[). Devel- 
oping microwave lattice clocks based on Cs may be ben- 
eficial as the active chamber of the clock will be reduced 
to a few micrometers across. This million-fold reduction 
in size is anticipated to lead to a better control over detri- 
mental black body radiation and stray magnetic fields. 

From Fig. [5] we find that the doubly-magic trap- 
ping of Cs atoms is indeed possible for two transitions: 
|4, 2) -> |3, -2) and |4, 3} -> |3, -3). The only complica- 
tion is that driving the former transition requires 4 pho- 
tons, while the latter transition requires 6 photons. This 
may be potentially accomplished either with multi-step 
RF/MW or stimulated Raman drives [9I [lrjj. 

The curve Q exhibits a resonant behavior when the 
laser frequency passes through the fine-structure doublet 
of atomic transitions: 6sx/2 — ®Pi/2 an d 6sx/2 — 6p 3 / 2 
at u = 0.050932 a.u. and 0.053456 a.u.. Not shown 
in the Fig. [2] are the values of the r.h.s. of Eq. © be- 
tween the resonances. In this region, the r.h.s. values be- 
come positive and the curve crosses the Mp>A cos 6 = +3 
limit from above: the 6-photon transition may be made 
doubly-magic when the laser is tuned to inside the fine- 
structure doublet. 

87 Rb serves as the secondary frequency standard. An 
analysis of magic conditions for this atom is carried out 
in Fig. [31 There is a single, \F' = 2, 1) ->• \F = 1, -1), 
transition of interest here. This is a two-photon transi- 
tion. Curiously, the Mp'AcosO curve nearly touches its 
limiting value at w w 0.0565 a.u. (A m ~ 806 nm) some- 
what below the 5si/2 — 5p 3 / 2 resonance. Here the r.h.s. 
of Eq. ([3]) reaches values of w —1.05, i.e., it is just 5% off 
the limiting value of -1. While not quite achieving the 
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FIG. 3: (Color online) Same as in Fig. [2] for 87 Rb transition 
\F' = 2,1) ->• \F = 1,-1) and for 85 Rb transitions \F' = 
3, 1) ->• \F = 2, -1) and \F' = 3, 2) ->• \F = 2, -2). 



To conclude, by working at the doubly-magic (Stark- 
Zeeman) conditions, one can greatly reduce sensitivity 
to spatial inhomogeneity due to trapping/bias fields and 
also reduce sensitivity to temporal fluctuations of the 
fields. It is anticipated that a variety of applications 
could take advantage of the magic (and nearly-magic) 
conditions found in this paper. For example, we antici- 
pate that lifetimes of quantum memory [171 ] may be im- 
proved. Another opportunity is developing microMagic 
lattice clocks. The important finding here is that multi- 
photon transitions in metrologically-important 133 Cs can 
be made simultaneously insensitive to both intensities of 
trapping lasers and also to fluctuations of magnetic fields. 
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"doubly-magic" status, this 806 nm wavelength gets us 
to nearly-magic conditions. In practice, if one may af- 
ford small Zeeman decoherences, the bias B-ficld may be 
detuned off its magic value. Then (see discussion above) 
the last term in the Stark shift, Eq.((5J|, is rescaled by the 
ratio (B/B m ) and the magic condition for the Stark shift 
can be reached. For example, I computed that A m ss 806 
nm becomes "Stark-magic" at B-field of 3.45 G, i.e., just 
6% larger than its "Zeeman-magic" value. Ultimately, 
the choice of B should be a matter of optimizing toler- 
ances to both Stark and Zeeman-induced decoherence in 
a particular application. For example, in a recent Paris 
experiment [13 ] with a magnetically trapped ensemble of 
87 Rb atoms, the bias B-field has been varied by as much 
as 16% from its magic value; this still has led to well- 
resolved contrast on the \F' = 2,1} — > \F = 1,-1) clock 
transition. 

Finally, in the same Fig.[3]we explore magic conditions 
for another isotope of Rb, 85 Rb. Compared to 87 Rb, the 
nuclear spin of this isotope is larger (J = 5/2), this re- 
sults in larger values of F and richer magnetic substruc- 
ture of the hyperfine levels. Two transitions of inter- 
est become available: \F' = 3, 1) — > \F = 2,-1) and 
\F' = 3,2) — > \F = 2,-2). Moreover, various contri- 
butions to Eq. @ scale differently with the nuclear spin 
and we need to carry out separate calculation for each 
isotope. From Fig. [3l we find that doubly-magic condi- 
tions can be attained for the \F' = 3,2) — > \F = 2, —2) 
transition in 85 Rb. 

I also carried out calculations for other commonly-used 
alkalis. For 7 Li, 23 Na, and 39 K there are no doubly-magic 
(or "near-magic") points; all these isotopes have 7 = 3/2. 
For various Fr isotopes, there is a multitude of doubly- 
magic points for multi-photon transitions. For example, 
for 210 Fr (I = 6) transition M F = -7/2 +7/2 these 
conditions are attained in the range A m = 846 — 1061 nm. 
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